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Abstract 

We compare the quantitative microstructural properties of Berea Sandstone with 
stochastic reconstructions of the same sandstone. The comparison is based on lo- 
cal porosity theory. The reconstructions employ Fourier space filtering of Gaussian 
random fields and match the average porosity and two-point correlation function of 
the experimental model. Connectivity properties of the stochastic models differ sig- 
nificantly from the experimental model. Reconstruction models with different levels 
of coarse graining also show different average local connectivity. 



Recently a number of stochastic models have been proposed for reconstruc- 
tion of the microstructure of porous media(see [1,2] and references therein). To 
assess the quality of the reconstruction, it is neccessary to have quantitative 
methods of comparison for such microstructures. General geometric character- 
ization methods normally include porosities, specific surface areas and correla- 
tion functions [4] . Here we follow a more general quantitative characterization 
for stochastic microstructures which is based on local porosity theory (LPT) 
[3,4]. Our analysis allows to distinguish quantitatively between three different 
microstructures all of which have identical porosities and correlation functions. 
The three microstructures are an experimental sample of Berea Sandstone ob- 
tained by computerized microtomography and two stochastic models of the 
same sandstone obtained through the Gaussian filtering method [1]. 

Consider a three-dimensional sample § = P U M (with P fl M = 0) where P is 
the pore space, M is the rock or mineral matrix. denotes the empty set. The 
porosity 0(§) of such a two component porous medium is defined as the ratio 
0(§) = V{¥) /V{^) where V{¥) denotes the volume of the pore space, and 
V(S) is the total sample volume. For the sample data analysed here the set § 
is a cube with sidelength M in units of the lattice constant a of a simple cubic 
lattice. Let K(r, L) denote a cube of sidelength L centered at the lattice vector 
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r. The set K(r, L) defines a measurement cell inside of which local geometric 
properties such as porosity and pore space connectivity are measured [3-5]. 
The local porosity in this measurement cell K(r, L) is defined as 0(r, L) — 
[V^(PnK(r,L))]/[y(K(r,L))]. The local porosity distribution n{(j),L) is defined 
as L) — — Y,r ^{4> ~ L)), where m is the number of placements of the 
measurement cell ]K(r, L). For better statistics the results presented here are 
obtained by placing ]K(r, L) on all lattice sites r which are at least a distance 
L/2 from the boundary of S. The local percolation probabilities characterize 
the connectivity of measurement cells of a given local porosity. Let A„(r, L) 
equal 1 if ]K(r, L) percolates in "a" direction and otherwise, be an indicator 
for percolation. A cell K(r, L) is called "percolating in the x-direction" if 
there exists a path inside the set P fl K(r, L) connecting those two faces of 
S that are vertical to the a;-axis. Similarly for the other directions. A3 = 1 
indicates that the cell can be traversed along all 3 directions, while Ac = 1 
indicates that there exists at least one direction along which the block is 
percolating. The local percolation probability in the "a" -direction is defined 
through Aa(0, L) = [Er A«(r, L)(500(r,L)]/[Er (^Mr.L)] and gives the fraction 
of measurement cells of size L having porosity that are percolating in the 
"a" -direction. The total fraction of percolating ce//s which percolate along the 
"a" -direction is given by Pa{L) — Jq jji{(f), L)Ao,(0, L) dcf). 

The Gaussian field (GF) reconstruction model [1] generates a random pore 
space configuration with inputs from a given experimental sample. Given the 
reference correlation function GEx(r) and porosity 0(Sex) of the experimental 
sample, the three main steps of constructing the sample §gf with correlation 
function G'gf(i") = G'ex(i") and porosity 0(Sgf) = 0(Sex) are as follows: 

(1) A standard Gaussian field X{r) is generated which consists of statistically 
independent Gaussian random variables X e R at each lattice point r. 

(2) The field X{r) is first passed through a linear filter which produces a 
correlated Gausssian field Y{r) with zero mean and unit variance. 

(3) The correlated field Y{t) is then passed through a nonlinear discretization 
filter which produces the reconstructed sample Sgf- 

For the process described in step 2, we have followed an alternate and compu- 
tationally more efficient method proposed in Ref. [1] that uses Fourier Trans- 
forms. An effective reconstruction requires a large separation ^ex ^ M where 
M is the sidelength(in pixels) of the sample and ^ex is the correlation length 
of the experimental reference, defined as the length such that GEx{r) ~ for 
r > ^EX- Violation of this condition leads to inaccuracy in the implementation 
of step 2 of the reconstruction, which in turn leads to a discrepancy at small 
r between Ggf(^) and GEx{r). This problem can be overcome by choosing 
large M. However, in d — 3 very large M also demands prohibitively large 
memory. Apart from this, the reconstruction also depends crucially on two 
other parameters, a length M(. up to which the experimental correlation is 
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incorporated into the reconstructed sample, and n, an interval at which the 
G'ex('") is sampled. For better reconstruction G'ex(^c) needs to be negligibly 
small. Different values of n correspond to a change of length scale. The model 
BRl is constructed with n — 1 and BR2 with n — 2. Although both have the 
same sidelength, the effective size of BR2 is twice that of BRl because of this 
coarse graining procedure. 

In Fig.l the averaged correlation functions G{r) = (G(r, 0,0) + G(0,r, 0) + 
^(0,0, r))/3 for the three samples are plotted. The experimental sample is 
a Berea sandstone of porosity = 0.1775 [5]. The model BRl, with 

M = 128, n = 1 and = 32, shows descrepancies for small values of r. 
However, for BR2, with M = 128, n = 2 and = 32, the correlation function 
matches more closely to that of EX. The resolution a of the experimental 
sample EX is 10/xm. Hence the actual size of EX and BRl is 1280//m, whereas 
that of BR2 is 2560/im. The porosities match quite well for all the samples 
(H^BRi) = 0.1783 and 0(§Bij2) = 0.1776). 

The reconstructed models BRl and BR2 are isotropic and globally connected, 
i.e., the pore spaces percolate in all the three directions. The local poros- 
ity analysis results are plotted in Fig. 2 for the experimental sample EX, 
in Fig. 3 for the stochastic model BRl, and in Fig. 4 for BR2. Compar- 
ison of L) indicates that the stochastic models have nearly the same 
level of homogeneity with that of the experimental sample. The main differ- 
ences are found in Aq(0, L) of the stochastic models. They differ significantly 
from that of EX, and they also vary widely among themselves. The recon- 
structed models have lower average connectivity of pore spaces. We observe 
[Xa{4>,L)]Ex > [)^a{4>,L)]BRi > [K{4>, L)]br2- Thesc differences appear 
even more clearly in the plot of Pa{L){mset of Fig 2, Fig 3 and Fig 4). In 
the experimental model(Fig 2) nearly all the measurement cells of dimension 
larger than 400/xm are percolating(globally connected pore space) in all di- 
rections. Comparison of p^{L) of the three models(Fig. 5) shows the drastic 
loss of average connectivity of the reconstructed models. Fig. 5 shows that 
the GF reconstruction BRl with n — 1 has a lower connectivity than EX. 
In BRl nearly 60% of the measurement cells of size 400//m percolate in all 
directions. Coarse graining [1] to n = 2 further reduces the connectivity of 
pore space. In BR2 less than 30% of the measurement cells of size 400;Lim 
percolate in all directions. These results indicate that gaussian filtering recon- 
struction methods retain a similar degree of isotropy and homogeneity as the 
original sandstone but tend not to reproduce connectivity properties that are 
important for transport. 



References 



3 



[1] p. Adler, Porous Media. (Boston: Butterworth-Heinemann, 1992). 

[2] C. Yeong and S. Torquato, Reconstructing random media, Phys. Rev. E 57 
(1998) 495-506. 

[3] R. Hilfer, "Local Porosity Theory for Flow in Porous Media", Phys.Rev. B 45 
(1992) 7115-7121. 

[4] R. Hilfer, Transport and relaxation phenomena in porous media. Advances in 
Chemical Physics XCII (1996) 299-424. 

[5] B. Biswal, C. Manwart and R Hilfer, Three-dimensional local porosity analysis 
of porous media, Physica A 255 (1998) 221-241 



4 




-0.2 I ' ' ' ' ' ^ 

50 100 150 200 250 300 

Fig. 1. Averaged directional correlation functions of all three models 




Fig. 2. Xa{(j),L) (broken curves, left axis) and fi ((f), L) (solid curve, right axis) at 
L = 200/xm for the model EX. The inset shows the function Pa{L), a is given in 
the legend. 
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Fig. 3. XaicpjL) (broken curves, left axis) and L) (solid curve, right axis) at 
L — 200/im for the model BRl. The inset shows the functions Pa{L), a is given in 
the legend. 




4> 

Fig. 4. Xa{<p,L) (broken curves, left axis) and /x(^,L) (solid curve, right axis) at 
L = 200/xm for the model BR2. The inset shows the functions Pa((L), a is given in 
the legend. 
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